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The  main  issue  in  preparing  fuel  cell  systems  for  the  future  market  is  system  reliability  and  efficiency.  Apart 
from  successful  field  test  trials,  any  type  of  stationary,  in  general  automotive  or  portable  fuel  cell  systems 
are  at  the  development  stage.  One  task  to  deal  with  is  to  increase  the  component  and  system  efficiencies  by 
facilitating  the  system  construction  or  eliminating  parasitic  components. With  newly  established  effective 
standardised  system  and  component  tests,  linked  with  a  flexible  modelling  and  simulation  environment, 
the  development  process  and  the  determination  of  the  system  efficiencies  as  well  as  the  inaccessible 
system  values  can  be  accelerated.In  this  work  a  modular  model-aided  system  analysis  and  development 
environment  is  presented  which  has  been  evaluated  and  validated  at  the  IWE.  The  tool,  a  combination  of 
standardised  testing,  modelling  and  simulation,  has  been  applied  to  different  types  of  fuel  cell  systems 
showing  the  tool  flexibility,  modularity  and  accuracy.  In  the  presented  case  the  tool  was  used  for  system 
analysis  and  studies  on  efficiency  increase  of  a  complex  prototype  stationary  PEMFC  system. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Standardised  testing,  modelling  and  simulation  play  key  roles 
in  the  improvement  of  fuel  cell  systems  and  applications.  The  sole 
realisation  of  hardware  tests  for  system  analysis  and  improvement 
is  not  enough,  often  takes  too  long  or  is  even  impossible.  During 
system  testing  many  important  values  and  a  plurality  of  parame¬ 
ters  necessary  for  a  detailed  analysis  or  even  an  exact  comparison 
of  different  systems  are  not  accessible.  As  a  consequence,  standard¬ 
ised  testing  has  to  be  linked  with  modelling  and  simulation  (Fig.  1 ), 
aiming  for  additional  component  and  system  information  and  a 
reduction  of  measurements.  The  limited  number  of  standard  test 
routines  will  therefore  deliver  the  necessary  model  parameters. 

At  the  IWE  a  model-aided  testing  and  simulation  approach  has 
been  evaluated  [  1,2  ]  and  extended  to  fit  all  kinds  of  fuel  cell  systems. 
For  validation  purposes  and  to  show  the  tool  flexibility,  modularity 
and  accuracy  for  a  wide  range  of  applications  it  has  been  applied 
to  a  complex  prototype  stationary  a  2  kWel  and  4  kWthermal  PEMFC 
Combined  Heat  and  Power  (CHP)  system  EDISon  [1-3]  (Fig.  2),  a 
1  kWel  commercial  methanol-based  PEMFC  Auxiliary  Power  Unit 
(APU)  and  a  Active  stationary  1  kWel  SOFC  CHP  system  [4]. 

In  this  work  the  tool  was  used  for  a  more  detailed  analysis  of 
the  prototype  stationary  2  l<Wel  PEMFC  system  EDISon.  The  results 
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of  the  performed  studies  on  possible  efficiency  increase  by  system 
design  improvement  and  optimised  system  operation  are  shown. 

2.  Standardised  system  testing 

Aside  the  reduction  of  system  tests  it  is  a  main  task  to  stan¬ 
dardise  them  in  order  to  get  comparable  system  test  results  and  to 
determine  additional  interesting  system  parameters  for  the  system 
and  component  models.  Therefore  a  harmonised  testing  format  has 
been  developed.  During  the  European  Commission  funded  project 
FCTESTNET  [5]  test  module  programs  have  been  evaluated  based 
on  existing  IWE  testing  experience  with  stationary  systems  [3]  or 
APU  systems  as  well  as  international  codes  and  standards  (Fig.  3). 
These  test  modules  are  going  to  be  further  proofed  in  the  European 
Commission  funded  project  FCTESQA  [6]. 

3.  Flexible  modelling  and  simulation  environment 

For  system  development  and  design,  many  companies  and 
research  institutes  use  models  which  fit  their  in-house  develop¬ 
ment.  But  for  an  effective  system  improvement  models  are  required 
flexible  to  fit  with  a  wide  variety  of  system  types.  A  further  impor¬ 
tant  attribute  of  a  system  modelling  and  simulation  tool  is  its  ability 
to  perform  dynamic  and  transient  simulation  studies  by  providing 
additionally  an  easy  to  handle  functionality. 

At  the  Universitat  Karlsruhe  (TH),  Institut  fur  Werkstoffe  der 
Elektrotechnik  (IWE)  a  flexible  and  easy  to  handle  dynamic  model- 
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Nomenclature 

A  area(m2) 

c  species  concentration  (mol  nrr3 ) 

Cp  molar  heat  capacity  (J  mol-1  K-1 ) 

d  thickness  (m) 

D  diffusion  coefficient  (m2  s-1 ) 

A rG  Gibbs  free  reaction  energy  (kj  mol-1 ) 

ArG°  Gibbs  free  reaction  energy  at  standard  pressure 
(kj  mol-1 ) 

A rH  reaction  enthalpy  (kj  mol-1 ) 

ArH  reaction  heat  flow  (W) 

A HVc  vaporisation,  condensation  reaction  heat  flow  (W) 

A rS  reaction  entropy  (kj  mol-1  K_1 ) 

Fact  activation  energy  (J,  eV) 

h  molar  enthalpy  (J  mol-1 ) 

H  enthalpy  (J) 

H  enthalpy  flow  (W) 

/  current  (A) 

j  current  density  (A  cm-2 ) 

jo  exchange  current  density  (A  cm-2 ) 

jFc  fuel  cell  current  density  (A  cm-2 ) 

j b  maximum  boundary  current  density  (A  cm-2 ) 

k  heat  transfer  coefficient  (W  nrr2  K-1 ) 

kA  area-averaged  heat  transfer  coefficient  (W  K-1 ) 

Kp  equilibrium  constant 

n  molar  amount  (mol) 

h  molar  flow  (mol  s_1) 

p°  standard  pressure  (1.013  bar) 

Pi  partial  pressure  of  the  component  i  (bar) 

psat  saturation  vapor  pressure  (bar) 

Q  heat  flow  (W) 

R  ohmic  resistance  (£2) 

R  gas  constant  (8.3145  J  mol-1 1<-1 ) 

Sco  percentage  CO  selectivity  (%) 

T  temperature  (K  or  °C) 

T0  temperature  at  standard  conditions  (298  K) 

Uc o  percentage  CO  conversion  (%) 

U  voltage  (V) 

U0  open  circuit  voltage  (V) 

Urev  reversible  cell  voltage  (V) 

x  extend  of  reaction  (mol  s-1 ) 

Greek  letters 

a  heat  conductivity  coefficient  (W  nr2 1<-1 ) 

8  thickness  (m) 

Pact  activation  overvoltage  (V) 

Pdiff  diffusion  overvoltage  (V) 

Pq  ohmic  overvoltage  (V) 

A  specific  heat  conduction  capability  (W  m_  1 1<-1 ) 

v  stoichiometric  coefficient 

p  specific  resistance  (Q  m) 


aided  analysis  environment  has  been  developed.  This  tool  is  able 
to  assist  hardware  test  routines  and  -  on  the  basis  of  the  accessible 
data  of  a  system  with  the  help  of  a  general  system  model  -  a  method 
has  been  validated  which  enables  an  evaluation  of  the  systems, 
helps  recognise  the  efficiency  of  individual  system  components, 
reveals  the  problems  during  operation  and  allows  an  appraisal  of 
the  further  optimisation  potential  of  the  systems. 

As  described  earlier  [1,2]  the  modelling  is  based  on  a  general 
variable  and  flexible  system  model  structure  consisting  of  dif¬ 
ferent  subsystem  groups.  This  structure  provides  a  frame  for  all 


kinds  of  fuel  cell  systems  (e.g.  combined  heat  and  power  systems 
(CHP)  or  auxiliary  power  units  (APU))  and  was  implemented  in 
Matlab/Simulink™,  creating  a  flexible  modelling  and  simulation 
environment. 


3.1  System  model  structure  and  topology 

The  main  benefit  of  the  presented  model-aided  analysis  tool  is 
its  modularity  and  standardised  build-up.  Starting  from  the  basis 
model  structure  the  functional  subgroups  are  split  into  further 
model  levels  according  to  the  chosen  top-down-strategy.  The  flexi¬ 
ble  model  structure  for  components  and  the  system  allows  the  use 
of  already  existing  and  own  developed  models  in  nearly  any  grade 
of  detail. 

In  Fig.  4  an  example  for  a  main  subsystem  group  (fuel  process¬ 
ing)  is  shown  including  the  standard  interface  arrays,  which  are 
used  in  the  subsystem  groups  and  the  models  itself  for  a  defined 
data  exchange. 

The  main  subsystem  groups  are  then  filled  by  connection  of 
different  component  models  matching  the  system  setup.  The  fuel 
processing  is  thus  able  to  represent  the  different  fuel  process¬ 
ing  setups,  e.g.  for  methane-based  stationary  PEMFC  systems, 
methanol-based  PEMFC  APU  or  methane-  and  diesel-based  station¬ 
ary  SOFC  systems  (Fig.  5). 

3.2.  Component  model  library  and  assumptions 

For  the  different  system  components  a  model  library  was 
established  (Fig.  6).  This  model  library  consists  of  steady  state 
or  dynamic  electrical  component  models  (blowers,  compressors, 
pumps,  inverters  and  other  power  electronic  devices),  general  com¬ 
ponent  models  including  a  species  and  heat  balance  for  steady 
state  and  dynamic  simulation  (reformer  steps,  heat  exchanger, 
vaporisers  and  burners)  and  extended  general  models  including 
electrochemistry  for  fuel  cells.  The  component  models  provide 
standardised  interfaces  and  allow  a  steady  state  and  dynamic  sim¬ 
ulation  of  systems  in  any  detail. 

The  standardised  model  structure  as  well  as  the  C-Code  com¬ 
patibility  of  Matlab/Simulink™  provides  an  easy  adaptation  of 
literature  models  [7,8]  or  models  in  more  detail  [9]  for  simulation. 


3. 2 A.  Chemical  equilibrium  calculation  and  species  balance 

For  the  chemical  equilibrium  calculation  and  the  species  balance 
ideal  gases,  incompressible  liquids,  total  mixing  and  no  pressure 
changes  in  the  reaction  volume  are  assumed. 

For  reversible  chemical  reactions  the  relation  between  the 
participating  species  and  their  stoichiometric  coefficients  v  is  as 
follows  [10,11]: 


uaA  +  vhB  ^  vcC  +  vdD 


(1) 


For  a  reaction  in  chemical  equilibrium  (ArG  =  0)  together  with 
the  definition  for  the  equilibrium  constant  I<p  it  can  be  written: 


In  I<p  = 


-A  rG°(T) 
RT 


(2) 


Eq.  (1 )  together  with  the  ideal  gases  law  leads  to  Guldberg’s  law 

[10,12]: 


„  h(C)v‘h(D)v«  1 


(3) 


The  extend  of  reaction  x  describes  the  converted  species  amount 
per  time  of  the  reactions  i.  With  the  stoichiometric  coefficients  of 
the  species  taking  part  at  the  reactions  and  the  conducted  molar 
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Model-Aided  System  Analysis  and  Development  Environment 


Standardized  System  Testing 


Physical  Modeling  and  Simulation 


Fig.  1.  Realised  solution  for  an  effective  analysis  and  optimisation  of  fuel  cell  systems. 


flows  hin  a  common  equation  system  according  to  Eq.  (3)  can  be 
established: 

/-A  rC° 

Kp,;  =  exP  y—Rj1- 

=  (n(Cj)in  +  VcjXjYfK  ri(Dj)in  +  Vd,iXj)Vd-t _ 1 _ 

(n(A-) in  +  Va,iXi)1'a-i(n(Bi)in  +  VbjXiY^  ^fl0^vc,i+vd,i-va,i~vb,i') 

(4) 

The  resulting  equilibrium  molar  flows  of  the  products  can  then 
be  written  as 


The  heat  of  the  reaction  flow  A rH  (exothermic  or  endothermic) 
can  be  calculated  out  of  the  sum  of  the  reaction  enthalpy  weighted 
with  extend  of  reactions: 

A  rH  =  Vx,ArHj  (6) 

i 

At  appropriate  pressure  and  temperature  conditions  vapori¬ 
sation  or  condensation  can  occur  in  the  reaction  volumes.  The 
resulting  heat  flow  A HVc  can  be  calculated  with  the  enthalpy  dif¬ 
ference  of  the  species  i  phase  change  from  liquid  (/)  to  gaseous  (g) 
and  vice  versa: 


"(AW  =  "(A)in  + 
"(AW  =  "(A) in  +  vb,i*i 
"(Q)  out  =  n(Ci)in  +  vc,i*i 
n(Dt)  out  =  ii(D,)in  + 


AHVC,i  =  Ahvc^r  =  (ftf  -  (7) 

(5)  The  enthalpy  difference  of  the  different  species  is  approximated 

by  temperature-dependant  polynomials  using  according  state  dia¬ 
grams  [13,14].  For  the  enthalpy  difference  of  water  the  following 


Fig.  2.  System  diagram  of  the  stationary  natural  gas  2  l<Wei  PEMFC  system  EDISon  at  the  IWE. 


International  Codes  &  Standards 


Harmonized  Testing  Format 


ASME  PTC  50-2002 

Fuel  Cell  Power  Systems  Performance 
ANSI/CSA  FC  1-2004 

Stationary  Fuel  Cell  Power  Systems  _ 

I  EC  62282-3-2 

Stationary  Fuel  Cell  Power  Plants  -  Test  Methods  for  the 
Performance 

DVGWVP  119 

German  code  for  Fuel  Cells  in  Combined  Heat  and  Power 
Generation  up  to  70  kW 


IWE  Testing  Experience 

i 

European  Projects: 
FCTESTNET 
FCTESQA 

i 

Test  Module  Program 


Fig.  3.  Definition  of  standardised  system  test  routines. 
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Fuel 


Air 


Cooling 

Process  Water 


Reformat-Offgas 

Electrical  Power 


Thermal  Power 


Fuel 

Processing 


Cooling 

Reformat 


Thermal  Power 


Waste  Water 


Exhaust 


-< 


Standard  interface  arrays 

•  np  total  mole  flow  (mol/s) 

•  p  pressure  (bar) 

•  T  temperature  (K) 

•  x_H2  mole  fraction  gas  or  fluid 

•  x_02 

•  x_CH4 

•  x_CO 

•  x_C02 

•  x_N2 

•  x_H20 

•  x_CH3OH 

•  x_CJ-i5OH 

•  x_gas/fluid_1 

•  P_therm  thermal  power  (W) 

•  P_el  electrical  power  (W) 

•  Qp  heat  flux  (W) 


Fig.  4.  Example  of  a  main  subsystem  group  including  the  composition  of  the  standardised  interface  array. 


polynomial  is  used  (T  in  °C): 

hH2o  -  hH2o  =  45070  -  41  -9T  +  3.44  x  10-3T2 

+  2.54  x  10“6T3  -  8.98  x  10“10T4  (J mol-1 )  (8) 

The  amount  of  the  vaporising  or  condensing  species  hf  is 
dependant  on  the  steam  mole  flow  hf  in  the  gas,  the  liquid  phase 
hi  and  the  temperature-dependant  saturation  amount  hf.  For 
the  heat  balance  modelling  a  prefix-afflicted  assignment  is  posted 
dependant  on  the  saturation  vapor  pressure  pf  and  the  total  pres¬ 
sure  ptot: 


3.2.2.  Energy  balance  and  heat  transfer 

For  the  energy  balance  and  the  heat  transfer  ideal  gases,  total 
mixing  and  incompressible  liquids  are  assumed. 

Volumes  and  masses  of  the  system  components  are  responsi¬ 
ble  for  the  dynamic  of  the  heat  transfer  influencing  the  thermal 
behaviour  of  the  components.  Regarding  energy  and  mass  conser¬ 
vation  for  any  open  system  a  dynamic  energy  balance  can  be  posted, 
which  includes  the  conducted  and  dissipated  heat  and  enthalpy 
flows  of  the  system: 

dH 

—  Qi  n  —  Qout  +  Ffin  —  f^out  (H  ) 


0, 


for{[(h 

for{[(h 


>  0)  A  (p 

>  0)  A  (p 


<  Ptot)]  A  [(hf  -  hf)  >  0]}  V  {[(hsat  >  0)  A  (Pf  <  ptot)]  A  [(-h\  <  hf  -  hf)  A  (hf  -  hf  <  0)]} 


<  Ptot)]  A  [(-h|  >  hf  -  hf)  A  (hf  -  hf  <  0)]}  V  {( pf  >  Ptot)} 


for  {[(hf  >  0)  A  (pf  <  Ptot)]  A  [(hf  =  0)  A  (hi  =  0)]} 


The  steam  saturation  amount  per  time  hf  of  a  species  in  the 
gas  flow  can  be  calculated  out  of  the  ideal  gases  law  and  the  partial 
pressures  and  molar  flows  of  the  different  components.  At  constant 
pressure  and  temperature  this  leads  to 


Pd 


d  H 


Ptot  -  P 


sat  "tot,  dry 


(10) 


dT 


v niCp,i  jyp  —  Qin  -  Qout  +  Hm  ~  ^out 


P’2  dt 


(9) 


According  to  [10]  the  enthalpy  H  of  ideal  gases  is 
H  =  nh  =  nCp(T-T°98K) 

Together  with  Eq.  (11 )  this  leads  to 


(12) 


(13) 


Fig.  5.  Different  types  of  fuel  processing  setups  represented  by  the  main  subsystem  group  fuel  processing. 
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Lookup  Table  (steady  state) 
and 

Electrical  Component  Models 


General 

Component 

Models 


Extended  General 
Component 
Models 


species  balance 
electrochemistry 
heat  balance 


Fig.  6.  Structure  of  the  component  model  library  providing  standardised  interfaces  for  steady  state  and  dynamic  simulations  of  systems. 


H[n  and  Hout  represent  the  sum  of  all  conducted  and  dissipated 
molar  enthalpy  flows  [10]: 


Hin  =  in  =  J>Cp,(T  -  r2°98K) 


1=1 


1  =  1 


Wout  —  ^  ^i, out  —  ^  ^jCp,i(T  ^298 1<) 
i=l  i=l 

Accordingly  for  the  heat  flows  Qin  and  QoUt: 


Qin  —  ^ 

i=l 

m 

Qout  =  ^  ^Qj, 


i,out 


i=l 


(14) 


(15) 


For  the  heat  conduction  through  a  thin  barrier  the  heat  flow 
is  dependant  on  heat  conduction  capability  X  of  the  material,  the 
thickness  d  of  the  barrier,  the  passed  through  barrier  area  A  and  the 
temperature  difference  AT  between  outside  and  inside  the  barrier: 


hydrogen  by  a  steam  reforming  unit  including  a  CO  shift  step.  A  CO 
purification  is  carried  out  by  selective  oxidation.  The  product  gas, 
which  is  fed  to  the  stack,  exhibits  a  CO  concentration  <30  ppm.  The 
PEMFC  stack  consists  of  80  cells  with  an  electrode  area  of  126  cm2 
each.  Further  components  in  the  system  are  a  methane  burner,  a 
desulphurisation  unit,  a  water  purifying  unit,  compressors  for  air 
and  fuel  and  a  power  inverter  for  grid  connection.  The  system  shall 
supply  max.  1.8  kW  electricity  and  approximately  4kW  heating 
power  for  domestic  hot  water  and  space  heating. 

The  system  was  used  for  performing  specific  simulation  studies 
regarding  hardware  and  efficiency  optimisation.  In  the  following 
chapters  the  main  system  components  model  assumptions  are 
mentioned. 

4.1.  Methane  steam  reformer  model  assumptions 

Together  with  the  steam  reforming  reaction  of  methane  the  side 
reaction  of  the  CO-shift  occurs: 

(1)  CH4  +  H2O^CO  +  3H2  ArH°9gK  =  206.1  kjmor1 

(2)  CO  +  H20  ^  C02  +  H2  ArH°9gK  =  -41.2  kjmol”1  1  J 


Q=^A7  (16) 

In  case  of  convection  a  gas  or  a  liquid  is  streaming  along  on 
side  of  the  barrier  absorbing  heat  and  dissipating  it  at  a  different 
place.  Analogue  to  Eq.  (16)  the  heat  flow  is  calculated  using  heat 
conductivity  coefficient  a: 

CL  =  aA  AT  (17) 


For  the  total  heat  transfer  through  a  barrier,  where  two  fluids 
(fluid  1  and  fluid  2)  are  streaming  along  (e.g.  a  heat  exchanger), 
the  convection  on  both  sides  as  well  as  the  heat  flow  through  the 
barrier  has  to  be  considered.  Therefore  a  heat  transfer  coefficient  k 
is  calculated  [15]: 


1  1  d  1 

T  ~ - F  y  “1 - 

k  ot\  X  a2 


(18) 


This  leads  to  the  heat  flow  transferred  from  volume  1  to  volume 

2: 


This  leads  together  with  Eq.  (4)  to  the  reaction  equation  system: 


Kp,i  =  exp 


-A rG°\  _  (hCojn+Xi  -X2)(hH2,in  +  3Xi  ±x2)3 


RT 


1 


(^tot,in) 


I<P,2  =  exp 


-A  rC° 
RT 


(«C H4,in  -  *1  )(nH20,in  “  *1  “  *2) 


±(21) 


(nco2,in  +  x2XnH2,in  +  3xi  +x2) 
(ncO,in+*l  -*2)(nH20,in-*l  ~  X2) 


The  resulting  equilibrium  molar  flows  of  the  product  gas  can  be 
calculated  out  of  Eq.  (5)  to 


nH2,out  =  nH2,in  +  3xi  ±x2 
HcO,out  =  nCo,in  +  x\  ~  x2 

t*C02,out  —  ^C02,in  +  *2  (22) 

HCH4,out  =  t*CH4,in  ~X\ 

^H20,out  =  t*H20,in  -  X\  -  X2 


Q  =  kA  AT  (19) 

The  heat  transfer  coefficient  and  the  temperature  difference  are 
only  valid  locally.  Therefore  area-averaged  heat  transfer  coefficients 
(/<A-values)  and  averaged  temperature  differences  (AT)  are  used  for 
calculation. 

4.  Stationary  PEMFC  CHP  system  process  model  and 
component  assumptions 

The  extensive  tested  and  analysed  stationary  PEMFC  CFIP  system 
EDISon  [1-3]  (Fig.  2)  operates  on  natural  gas,  which  is  converted  to 


The  simulation  parameters  (reaction  volumes  and  heat  transfer 
coefficients)  for  the  further  component  models  of  the  reforming 
unit  (steam  reformer,  burner,  vaporiser  and  heat  exchanger)  are 
listed  in  Table  1. 

4.2.  CO-shift  model  assumptions 

The  model  for  the  medium  temperature  shift  (MTS)  consists  of 
the  reaction  equation  for  the  CO  shift: 

CO  +  H20  ^  C02  +  H2  ArH°98K  = -41.2  kj mol"1  (23) 


126 


P.  Konig,  E.  Ivers-Tiffee  /  Journal  of  Power  Sources  190  (2009)  121-132 


Table  1 

EDISon  CHP  reaction  volumes  and  heat  transfer  coefficients— reforming  unit 


Component 

Volume  primary 

(m3) 

kA-value  primary 
(WK4) 

Volume 
secondary  (m3) 

kA-value 

secondary  (WK-1) 

Mass  solid  (kg) 

kA-value  solid 
(WK-1) 

Methane  steam  reformer 

4.37e-3 

7 

- 

- 

6 

0.19 

Methane  burner 

3.69e-3 

0.1294 

- 

- 

3 

0 

Vaporiser 

9.13e-4 

2.19 

le-4 

3.1 

0.703 

1.2 

Heat  exchanger 

2.54e-3 

1.2 

1.3e-4 

2.5 

0.824 

0.3 

Exhaust  heat  exchanger 

1.81e-4 

2.7 

6.41  e-3 

6.3 

1.144 

1.512 

The  resulting  equilibrium  molar  flows  of  the  product  gas  are 


The  resulting  equilibrium  molar  flows  of  the  product  gas  then 
are 


^H2,out  —  ^H2,in  +  * 


f*CO,out  =  f*CO,in  “  x 


f*C02,out  —  f*C02,in  T  x 


(24) 


f*H2,out  —  ^H2,in  2[ho2jn  ^02,out  0.5(hco,in  ^CO,out)] 

Uco  . 

nCO,out  =  nCO,in  -  ^QQnCO,in 


f*H20,out  —  ^H20,in  —  x 

The  simulation  parameters  (reaction  volumes  and  heat  transfer 
coefficients)  for  the  further  component  models  of  the  CO- 
shift  unit  (vaporiser,  CO-shift  and  heat  exchanger)  are  listed  in 
Table  2. 


4.3.  CO  purification  model  assumptions 


flC02,out  =  ^C02,in  +  (^CO,in  “  HcO.out)  (27) 

HH20,out  =  f*H20,in  +  2[ho2,in  “  ^02,out  -  0.5(hco?in  -  hco,out)] 
(^CO,in  “  nCO,out)100 

n02’Out  -  "°2.in  2Sco 

The  simulation  parameters  (reaction  volumes  and  heat  transfer 
coefficients)  for  the  further  component  models  of  the  CO  cleaning 
unit  (SelOx  and  heat  exchanger)  are  listed  in  Table  3. 


For  the  CO  purification  using  a  selective  methanation  (SelMet) 
the  model  for  the  methane  steam  reforming  can  be  used  as  reverse 
reaction. 

For  the  CO  purification  using  a  selective  oxidation  (SelOx)  aside 
the  CO  oxidation  the  H2  oxidation  occurs  in  parallel: 

(1) CO+  ^02  ^C02  ^rH298K  =  -282.9 kj mol-1 

(25) 

(2) H2  +  lo2  -  H20  ArH°98K  =  -241. 8  kj  mol-1 


4.4.  PEMFC  stack  model  assumptions 

For  the  fuel  cell  model  ideal  gases,  constant  pressures  along  the 
channels,  electrolyte  only  proton  conductive  and  no  further  fluid 
transport  are  assumed. 

The  reaction  equation  for  the  H2-oxidation  is 

H2  +  io2^H20  ArH^<  = -285.8  kjmor1 2  (28) 


Setting  a  percentage  CO  conversion  Uc o  and  a  percentage  CO 
selectivity  Sco  the  model  can  be  adapted  to  different  real  reformers 
and  measurement  values: 


The  hydrogen  requirement  hH2,req  of  the  fuel  cell  can  be  calcu¬ 
lated  depending  on  load  current  and  number  of  cells  using  Faraday’s 
law  [16]: 


UC  o  = 


Sco  = 


nCO,  utilised 
^CO,in 

^02 ,  CO-oxidation 
^02,  utilised 


(26) 


N  NAce ii  . 

nH2,req  =  ^load  =  (29) 

The  output  molar  flow  composition  of  the  anode  and  the  cath¬ 
ode  are  calculated  using  Eq.  (28)  and  the  chemical  equilibrium 


Table  2 

EDISon  CHP  reaction  volumes  and  heat  transfer  coefficients— CO-shift  unit 


Component 

Volume  primary 

(m3) 

kA-value  primary 
(WK-1) 

Volume 
secondary  (m3) 

kA-value 

secondary  (WK-1) 

Mass  solid  (kg) 

kA-value  solid 
(WK-1) 

Vaporiser 

0.05 

- 

- 

- 

- 

- 

CO-shift 

0.88e-6 

1.25 

- 

- 

0.5 

0.12 

Heat  exchanger 

0.52e-3 

6.5 

0.52e-3 

19 

le-5 

6 

Table  3 

EDISon  CHP  reaction  volumes  and  heat  transfer  coefficients— CO  cleaning  unit 


Component 

Volume  primary 

(m3) 

kA-value  primary 
(WK-1) 

Volume  secondary 

(m3) 

kA-value 

secondary(WI<-1) 

Mass  solid  (kg) 

kA-value  solid 
(WK-1) 

Gas  to  gas  heat  exchanger 

1.32e-3 

2 

5.21e-3 

2 

le-5 

2 

Gas  to  liquid  heat  exchanger 

3.77e-6 

3 

3.39e-3 

2.8 

0.15 

0.01 

SelOx 

le-5 

14 

- 

- 

le-5 

1 

P.  Konig,  E.  Ivers-Tiffee  /  Journal  of  Power  Sources  190  (2009)  121-132 


127 


calculation  regarding  vaporisation  and  condensation: 


t*H2Aout  —  ^H2,A,in  —  ^H2,req 

ng  —  ng  —  nvc 

aH20,A,out  ~~  aH20,A,in  aH20 ,A 

1  (30) 
^02,C,out  —  t*02,C,in  —  2^H2’req 

^H20,C,out  =  ^H20,C,in  +  ^H2,req  “  ^ 0,C 

The  cell  voltage  is  assumed  to  be 

t/cell  =  ^rev  -  ^7act  -  ^diff  “  r 1 £2  (31 ) 


The  temperature  and  pressure  dependant  reversible  cell  voltage 
for  the  hydrogen  oxidation  can  be  written  as  [10,16-19]: 


Urev  — 


-A  rG° 
zF 


/  PH2y/PO^\ 

\Ph2oVp°  J 


(32) 


The  exchange  current  density  j0  is  temperature  and  concen¬ 
tration  dependant  [16,17].  In  [19]  the  exchange  current  density  is 
described  dependant  on  the  concentrations  Cox  and  Cred,  the  activa¬ 
tion  energy  Eact,  the  temperature  T  and  a  so-called  frequency  factor 
A: 


j0  =  zFAc'0-ac"ed  e-W«-  (33) 

In  the  model  the  activation  overvoltage  rjac t  is  calculated  regard¬ 
ing  an  anode  and  cathode  part  for  the  exchange  current  density  in 
Eq.  (33).  Assuming  constant  concentrations  at  anode  and  cathode, 
these  concentrations  lead  together  with  the  frequency  factor  to  a 
new  factor  A.  The  activation  energy  Eact  and  the  Boltzmann  constant 
are  combined  to  a  factor  B.  In  case  of  ideal  gases  the  concentrations 
can  be  replaced  by  partial  pressures  pz.  This  leads  to  the  expres¬ 
sions  for  the  exchange  current  densities  at  the  anode  j0  A  and  at  the 

cathode  jo, c- 


Jo,A=A^““A)|^e  Ba/Ta 

Jo,c-Acp^§-e-Bc/rc 


(34) 


For  the  modelling  of  the  diffusion  a  simplified  one-dimensional 
system  assuming  a  linear  concentration  decrease  through  the  elec¬ 
trode  is  used  [16,19].  Together  with  the  first  Fields  law  it  can  be 
written  as 


1  ;FC  ,pn  Ci  ~  CPB,i 

JFC  =  7 -  =  ZtUi - X - 

Acell 


(35) 


The  current  limit  is  reached,  when  the  boundary  phase  reaction 
is  faster  than  the  charge  carrier  diffusion  through  the  electrodes 
and  thus  reaching  a  zero  phase  boundary  concentration  cPBiI.  This 
leads  to  the  maximum  boundary  current  density  jB.max  using  the 
ideal  gases  law: 


Js.max  =  zFDj  |  =  zFD,  (36) 

At  the  boundary  between  gas  volume  and  electrode  the  potential 
U is  [19]: 


Ur 


At  the  phase  boundary: 

jr  RT  cPBi 
PB,i  =  u0  -  tfln~jr 


(37) 


(38) 


Table  4 

PEMFC  stack  voltage  current  curve  model  parameters 


Open  circuit  voltage 

U0 

0.95  V 

Fuel  utilisation  anode 

- 

90% 

Air  utilisation  cathode 

- 

30% 

Exchange  current  density  anode 

JO, A 

280  mA  cm-2 

Exchange  current  density  cathode 

J0,K 

1.8  mA  cm“2 

Transfer  factor  anode 

aA 

0.4 

Transfer  factor  cathode 

UK 

0.72 

Specific  resistance  anode 

Pa 

1.4e-3  Q  cm 

Specific  resistance  membrane 

Pm 

25  £2  cm 

Specific  resistance  cathode 

Pk 

1.2e-3£2cm 

Thickness  anode 

5a 

100  |jim 

Thickness  membrane 

5m 

175  |Jim 

Thickness  cathode 

5k 

100  fxm 

Diffusion  coefficient  anode 

Da 

le-7  m2  s_1 

Diffusion  coefficient  cathode 

D  K 

9e-8m2  s_1 

Factor  anode 

Ba 

0 

Factor  cathode 

Bk 

1/5000 

Table  5 

EDISon  CHP  reaction  volumes  and  heat  transfer  coefficients- 

-PEMFC  stack  unit 

kA-value  anode 

7WIC1 

kA-value  cathode 

7WK-1 

kA-value  cooling 

50  WK4 

kA-value  solid 

0.5WK-1 

Mass  solid 

12  kg 

Volume  cooling 

0.51 

Channel  size  anode 

0.5  mm 

Channel  size  cathode 

0.5  mm 

The  difference  of  the  electrochemical  potentials  within  the  elec¬ 
trodes,  caused  by  the  species  transport,  leads  to  the  diffusion 
overvoltage  pdiff: 


RT  cPB  i 

Bditt  =  U  -  UPB  i  = 


=  ,-JgL 

zF  \  zFDiCj 


RT , 

=  —pin 
zF 


Jfc  A 

jB,max  J 


(39) 


For  the  Ohmic  overvoltage  tjq  a  linear  increase  with  the  current 
density  is  assumed.  Therefore  constant  specific  resistances  pz  for 
the  electrodes  and  for  the  membrane  can  be  implied: 


t\q  =  Rkc  =  (Pa^a  +  Pm^m  +  Pc^c)Jfc  (40) 

The  stack  voltage  current  curve  model  parameters  used  in  the 
simulations  are  displayed  in  Table  4. 

Additional  simulation  parameters  (geometry,  reaction  volumes 
and  heat  transfer  coefficients)  for  the  PEMFC  stack  unit  are  listed 
in  Table  5. 


Table  6 

EDISon  CHP  steady  state  operating  points,  coolant  25  °C,  fuel  utilisation  50% 


System  operating  point 

50% 

60% 

70% 

Conducted  methane  heating  power  (kW) 

6.65 

8.41 

8.5 

Produced  H2  power  (kW) 

3.7 

4.4 

5.15 

Electrical  power  (kW) 

0.935 

1.09 

1.15 

Thermal  power  (kW) 

1.5 

1.9 

2.15 

Electrical  efficiency  (%) 

14.06 

12.96 

13.53 

Electrical  net  efficiency  (%) 

8.05 

8.2 

8.82 

Thermal  efficiency  (%) 

22.6 

22.59 

25.3 

Total  net  efficiency  (%) 

30.65 

30.79 

34.12 

Total  efficiency  (%) 

36.66 

35.55 

38.83 
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Fig.  7.  Response  dynamics  for  measured  and  simulated  electrical  system  power  output. 


Fig.  8.  Response  dynamics  for  measured  and  simulated  thermal  system  power  output. 


5.  Results  and  discussion 

The  measurement  results  of  the  1.8  kWel  and  4  kWthermal  PEMFC 
CHP  prototype  system  EDISon  showed  a  limited  possible  system 
operating  range  of  50-70%  maximum  electrical  output  power  due 
to  the  system  configuration  and  control.  The  suboptimal  cathode 
flow  duct  concept  of  the  PEMFC  stack  leads  to  a  maximum  possi¬ 
ble  fuel  utilisation  of  only  50%  over  the  system  operating  range.  In 
Table  6  the  reached  steady  state  system  operating  points  at  25  °C 
coolant  temperature  and  at  a  fuel  utilisation  of  50%  are  displayed. 

Using  the  developed  physical  models  for  the  different  reformer 
steps,  PEMFC,  heat  exchangers  and  compressors  a  system  model 
was  built  to  simulate  the  different  testing  procedures  performed 
prior  to  the  simulation.  In  Figs.  7  and  8  a  measured  typical  stan¬ 
dardised  dynamic  test  module  is  shown.  During  the  measurement 
the  system  output  power  is  varied  by  changing  the  load  current 
including  an  inverter  shutdown.  The  simulation  results  show  a  good 
agreement  with  the  measured  electrical  and  thermal  power  out¬ 


puts.  Even  the  electrical  load  step  response  dynamics  of  the  PEMFC 
stack  power  output  can  be  simulated  with  a  good  accuracy  (Fig.  7). 

Based  on  the  measurement  results  the  former  validated  model- 
aided  system  analysis  environment  is  now  used  to  detect  the 
optimisation  potential  of  the  prototype  PEMFC  CHP  system  and  to 
show  the  available  real  performance  values  of  the  system  and  its 
components. 

The  identified  main  criteria  resulting  in  poor  system  per¬ 
formance  are:  power  loss  due  to  falsely  designed  and  faulty 
components,  efficiency  loss  due  to  parasitic  auxiliaries,  efficiency 
loss  due  to  ineffective  component  setup  and  system  technology, 
poor  system  insulation  and  heat  integration  as  well  as  inappropri¬ 
ate  concepts  for  fuel  treatment  and  fuel  cell  operation.  All  these 
points  are  directly  influencing  the  efficiency  and  therefore  the 
performance  of  the  systems.  Aside  the  economical  optimisation 
and  the  long-term  stability  the  increase  of  the  efficiency  there¬ 
fore  represents  a  main  development  demand  for  stationary  fuel 
cell  systems. 
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Fig.  9.  Simulation  results:  thermal,  electrical  and  total  efficiency  of  the  stationary  natural  gas  21<Wei  PEMFC  system  without  (light  colour)  and  with  (dark  colour)  anode 
offgas  recirculation  to  the  burner. 


5.1.  Increase  of  the  fuel  utilisation 

The  measurement  results  and  the  performed  post-test  simula¬ 
tion  analyses  on  the  PEMFC  CHP  system  showed,  that  for  a  as  high 
as  possible  efficiency  the  fuel  stack  as  the  central  energy  converter 
has  to  be  designed  and  operated  as  optimal  as  possible.  Especially 
during  reformate  gas  operation  the  maximum  achievable  fuel  util¬ 
isation  in  the  fuel  cell  stack  has  an  essential  influence  on  system 
performance  and  efficiency. 

The  maximum  achievable  system  efficiencies  at  varying  fuel 
utilisation  are  listed  in  Table  7.  By  increasing  the  fuel  utilisation 
from  50%  to  an  aiming  realistic  value  of  80%,  absolute  efficiency 
increases  of  5.4%  electrical  and  11.5%  thermal  have  been  deter¬ 
mined  for  the  analysed  system.  This  can  be  achieved  by  optimising 
the  flow  duct  in  the  stack  and  reducing  the  inert  gases  in  the 
reformat  gas  flow. 

5.2.  Reduced  electric  auxiliaries  consumption 

In  every  steady  state  operating  point  the  analysed  system  has 
an  electrical  power  consumption  due  to  the  auxiliaries  of  about 
400  W.  In  relation  to  the  produced  power  there  is  a  high  poten¬ 
tial  for  efficiency  optimisation.  In  Table  8  the  achievable  system 


efficiencies  at  a  fuel  utilisation  of  70%  and  an  averaged  auxiliaries 
power  consumption  of 400  W  are  displayed.  By  reducing  the  power 
consumption  down  to  100  W  an  electrical  net  efficiency  increase  up 
to  3.5%  can  be  achieved. 

5.3.  Reduced  primary  energy  consumption— feeding  back  of  the 
anode  offgas 

Analysing  the  performed  system  measurements  and  simula¬ 
tions  [  1,2  ]  it  was  possible  to  quantify  the  feasible  efficiency  increase 
of  the  methane-based  stationary  PEMFC  CHP  system  when  feeding 
back  the  anode  offgas  to  the  burner  (Fig.  2)  and  therefore  saving 
methane.  In  Fig.  9  the  simulation  results  predict  an  average  absolute 
efficiency  increase  of  5%  electrical  and  3%  thermal. 

5.4.  Reduced  primary  energy  consumption— increase  of  the 
reformer  heat  insulation 

Using  an  infrared  camera  the  radiated  heat  losses  of  the  PEMFC 
CHP  system  reformer  unit  have  been  visualised.  At  the  specific 
system  operating  point  (50%  system  power)  a  reformer  outlet  tem¬ 
perature  of  803.4  °C  and  a  methane  consumption  in  the  burner  of 
4.604  x  10-3  mol  s-1  have  been  measured.  The  average  heat  trans¬ 


late  7 

EDISon  CHP  efficiencies,  variation  of  fuel  utilisation,  coolant  25  °C,  power  setpoint 
70% 


Fuel  utilisation  anode 

50% 

70% 

80% 

Electrical  efficiency  (%) 

12.99 

16.84 

18.37 

Electrical  net  efficiency  (%) 

7.65 

11.51 

13.04 

Thermal  efficiency  (%) 

21.17 

29.21 

32.65 

Total  efficiency  (%) 

34.16 

46.05 

51.02 

Total  net  efficiency  (%) 

28.82 

40.72 

45.69 

Table  8 

EDISon  CHP  variation  of  electric  losses,  coolant  25  °C,  power  setpoint  70% 

Electrical  auxiliaries  losses 

400  W 

200  W 

100W 

Conducted  methane  heating  power  (kW) 

8.5 

8.5 

8.5 

Electrical  power  (kW) 

1.431 

1.431 

1.431 

Electrical  efficiency  (%) 

16.84 

16.84 

16.84 

Electrical  net  efficiency  (%) 

12.13 

14.48 

15.66 

4.7E-03 

4,6E-03 

4.5E-03 

4.4E-03 

4.3E-03 

4.2E-03 

4.1E-03 

4.0E-03 

3.9E-03 

3,8E-03 


Fig.  10.  Simulated  possible  methane  consumption  reduction  of  the  burner  by 
increasing  the  insulation  as  well  as  the  resulting  increase  of  the  total  system  effi¬ 
ciency  at  50%  system  power. 
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Table  9 

EDISon  CHP  simulated  gas  composition  after  the  CO-shift  step 


Parameter 

Gas  composition 

Educt  H20 

0.02976  mol  s”1 

H2  (vol.%) 

79.43 

Educt  CH4 

0.00473  mols-1 

CO  (vol.%) 

0.3063 

Educt  N2 

0.00116  mol  s_1 

C02  (vol.%) 

19.63 

Pressure 

1.04  bar(a) 

CH4  (vol.%) 

0.1255 

Temperature  reformer 

806.9  °C 

N2  (vol.%) 

0.5143 

Temperature  CO  shift 

275.6  °C 

fer  coefficient  (kA- value)  from  the  reformer  solid  mass  to  the 
environment  has  then  been  determined  by  simulation. 

The  model  was  then  parameterised  with  these  values.  With  sim¬ 
ulation  studies  the  possible  decrease  of  the  heat  losses  by  increasing 
the  reformer  insulation  was  determined.  Fig.  10  shows  the  reduc¬ 
tion  of  the  natural  gas  molar  flow  depending  on  the  insulation 
as  well  as  the  resulting  total  system  efficiency.  By  increasing  the 
insulation  the  necessary  natural  gas  molar  flow  is  reduced.  The 
amount  for  the  reforming  is  kept  constant.  Using  a  three  times 
larger  insulation  (kA- value:  0.063  W  K-1 )  the  necessary  natural  gas 
amount  would  be  reduced  by  11.6%.  This  would  then  result  in  an 
absolute  system  efficiency  increase  (excluding  auxiliaries  losses)  of 
2.45%. 

5.5.  Optimised  concept  for  the  fuel  treatment 

An  essential  advantage  of  the  developed  modelling  and  simula¬ 
tion  approach  is  the  possible  simulation  of  different  kinds  of  system 
configurations  by  easily  replacing  components  in  the  model. 

A  very  promising  point  in  system  optimisation  is  to  rearrange  or 
to  exchange  system  components  to  improve  the  thermal  integra¬ 
tion  or  to  reduce  parasitic  components  and  to  simplify  the  system 
technology.  The  present  natural  gas  reformer  unit  of  our  Station¬ 


Table  n 

EDISon  CHP  system  types  1  and  2,  determined  maximum  system  efficiencies 


System  type 

SelOx 

SelMet 

Electrical  efficiency  (%) 

16.84 

17.01 

Electrical  net  efficiency  (%) 

11.51 

11.67 

Thermal  efficiency  (%) 

29.21 

28.95 

Total  efficiency  (%) 

46.05 

45.96 

Total  net  efficiency  (%) 

40.72 

40.62 

ary  PEMFC  system  consists  of  a  subunit  including  a  steam  reformer 
step  (STR)  and  a  CO  cleaning  step  (MTS).  The  needed  CO  fine  clean¬ 
ing  step  (selective  oxidation)  is  installed  separately  including  the 
needed  additional  air  dosing  increasing  the  auxiliary  losses. 

To  be  able  to  remove  the  necessary  air  dosing  and  compressed 
air  supply,  in  our  model  the  SelOx  has  been  replaced  with  a  selective 
methanation  step  (SelMet)  (Fig.  11). 

Additionally  a  better  thermal  integration  can  be  achieved  by  a 
direct  mounting  of  the  SelMet  into  the  reformer  unit.  The  exother¬ 
mic  reaction  of  the  SelMet,  as  reverse  reaction  of  the  methane  steam 
reformation,  allows  an  optimised  thermal  management  and  there¬ 
fore  saving  reformer  unit  heating  power. 

The  main  important  aspect  was  to  guarantee  that  the  gas  quality 
of  the  resulting  reformation  was  equal  or  even  better  than  before. 
Performed  simulation  studies  were  able  to  proof  the  concept.  In 
Table  9  the  selected  model  parameters  of  the  steam  reformer  and 
the  CO  shift  as  well  as  the  resulting  output  gas  composition  of  the 
CO  shift  at  60%  system  power  are  displayed. 

In  Table  10  the  according  resulting  product  gas  compositions  for 
the  SelOx  and  SelMet  at  the  60%  system  operating  point  are  listed. 

For  the  SelMet  the  simulation  shows  a  higher  hydrogen  amount 
of  79.22%  in  the  product  gas  and  nearly  no  CO  compared  to  the 
SelOx.  Due  to  the  different  chemical  reactions  the  SelMet  shows 
a  lower  methane  conversion  rate.  In  comparison  to  the  reduced 


Table  10 

EDISon  CHP  gas  quality  comparison  SelOx  and  SelMet,  power  setpoint  60% 


Parameter 

SelOx 

SelMet 

Gas  composition 

SelOx 

SelMet 

Temperature 

45  °C 

246.85  °C 

H2  (vol.%) 

75.09 

79.22 

Pressure 

1.04  bar(a) 

1.04  bar(a) 

CO  (ppm) 

5.33 

2.364  x  10“5 

Educt  02 

0.00148  mols”1 

- 

C02  (vol.%) 

19.28 

19.5 

Utilisation  CO 

99.82% 

- 

CH4  (vol.%) 

0.1214 

0.7439 

Utilisation  C02 

- 

1.5% 

N2  (vol.%) 

5.248 

0.5269 

CO  selectivity 

14.5% 

- 

02  (vol.%) 

0.2591 

- 

CH4  conversion  (%) 

99.4 

96.38 
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Table  12 

Optimised  EDISon  CHP  system  type  2  compared  to  pre-commercial  systems 
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system  components  (no  additional  air  input)  and  the  optimised  CO 
reduction  yet  this  is  negligible. 

The  small  amount  of  nitrogen  in  the  product  gas  of  the  SelMet 
is  a  result  of  the  conducted  natural  gas  composition.  This  obvious 
minimised  inert  amount  of  nitrogen  in  the  conducted  reformate 
gas  flow  of  the  SelMet  in  comparison  to  the  SelOx  reduces  the 
occurring  species  transport  blockings  at  high  load  currents  in  the 
PEMFC  stack  and  therefore  better  hydrogen  conversion  rates  are 
possible. 

The  absence  of  the  SelOx  in  the  coolant  circuit  of  the  fuel  cell 
stack  unit  also  allows  higher  stack  operating  temperatures  and 
therefore  higher  power  levels  due  to  the  missing  limited  SelOx 
operating  temperature  of  45  °C. 

The  benefit  of  the  possible  efficiency  increase  is  only  then  valid, 
if  already  at  the  same  operating  conditions  (coolant  temperature, 
fuel  utilisation,  load  current)  as  with  the  SelMet  similar  output 
power  levels  are  reached.  The  simulation  results  of  the  power  out¬ 
put  comparison  of  the  different  system  configurations  with  SelOx 
(EDISon  CHP  system  type  1)  and  SelMet  (EDISon  CHP  system  type 
2)  are  displayed  in  Table  11. 

The  SelMet  system  achieves  a  higher  electrical  efficiency  of 
0.17%  absolute.  The  SelOx  system  has  therefore  a  0.26%  absolute 
higher  thermal  efficiency.  The  total  efficiencies  of  both  systems 
are  therefore  almost  identical  and  the  above-mentioned  advan¬ 
tages  of  the  EDISon  CHP  system  type  2,  better  thermal  integration 
and  absence  of  the  air  supply,  allow  an  effective  system  efficiency 
increase. 


5.6.  Efficiency  optimised  simulated  system  compared  to  real 
pre-commercial  systems 

Out  of  the  sum  of  the  above-mentioned  possible  efficiency 
increases,  determined  with  the  standardised  test  procedures  and 
the  modelling  and  simulation  approach,  a  minimum  value  pre¬ 
diction  for  the  achievable  system  efficiencies  can  be  stated  and 
compared  to  existing  systems  (Table  12). 

Therefore  the  EDISon  CHP  system  type  2  (SelMet)  serves  as  an 
optimisation  basis  regarding  the  efficiency  increase  using  anode 
recirculation  and  optimised  reformer  insulation.  Additionally  the 
possible  absolute  efficiency  increases  achieved  with  a  fuel  utilisa¬ 
tion  increase  from  70%  to  80%  (simulated  values:  1.53%  electrical, 
3.44%  thermal)  as  well  with  a  reduction  of  the  auxiliaries  electrical 
consumption  to  100  W  (3.53%  net  electrical,  Table  8)  are  considered. 

Comparing  the  predicted  values  of  the  EDISon  2  system  with 
the  data  of  the  almost  similar  concept  of  the  industrial  Viessmann 
lab  prototype  2  [20,21]  the  applicability  and  ability  for  realistic 
and  exact  predictions  of  the  developed  model-aided  analysis  and 
simulation  approach  can  be  proven. 

Having  a  still  optimised  reformer  concept  the  Viessmann  lab 
prototype  achieves  a  better  total  system  efficiency.  With  the  addi¬ 
tional  improvements  by  a  more  compact  reformer  concept,  lower 
reformer  temperatures  and  the  resulting  reduction  of  primary 
energy  similar  to  the  Viessmann  lab  prototype  3  [20,21  ]  the  EDISon 
2  system  could  achieve  comparable  electrical  system  efficiencies  of 
about  30%. 


6.  Conclusions 

A  model-aided  system  analysis  and  optimisation  environment 
based  on  Matlab/Simulink™  has  been  developed  and  validated  for 
different  kinds  of  fuel  cell  systems.  The  approach  of  the  standard¬ 
ised  modelling  and  simulation  environment  including  the  likewise 
developed  standardised  system  testing  procedures  [5,6]  provides 
additional  information  about  input  and  output  values,  internal  vari¬ 
ables  of  state  as  well  as  the  system  efficiency  in  the  steady  state  and 
in  the  dynamic  operation  mode  and  thus  proves  itself  a  powerful 
tool. 

The  knowledge  of  the  power  and  operating  limits  of  the  fuel  cell 
systems  leads  to  the  development  of  optimisation  strategies  and 
the  determination  of  the  realisable  system  and  component  poten¬ 
tial.  In  this  work  the  most  important  points  for  the  judgement  of 
the  system  and  component  power  and  optimisation  potential  have 
been  analysed: 

•  efficiency,  modulation  ability  and  system  operating  range; 

•  application  evaluation  of  alternative  or  additional  system  com¬ 
ponents; 

•  optimisation  of  the  thermal  integration  and  reduction  of  thermal 
losses; 

•  hardware  optimisation:  system  setup,  reduce  of  parasitic  system 
components. 

Aside  the  economical  optimisation  and  the  long-term  stability 
the  increase  of  the  system  efficiency  represents  the  main  develop¬ 
ment  demand  for  stationary  fuel  cell  systems.  Using  the  developed 
modelling  and  simulation  approach  as  well  as  the  developed 
standardised  test  procedures,  a  prototype  of  a  stationary  natural 
gas-based  2  kWel  PEMFC  CHP  system  has  been  analysed  regard¬ 
ing  power  and  optimisation  potential.  Thereby  the  most  important 
criteria  have  been  identified  influencing  the  PEMFC  CHP  system 
efficiency. 

Aside  a  preferably  high  fuel  utilisation  rate  the  reduction  of  the 
auxiliaries  electric  power  consumption  is  one  of  the  main  tasks. 
By  reducing  the  electrical  power  consumption  down  from  400  to 
100  W  an  electrical  net  efficiency  increase  up  to  3.5%  could  be 
achieved. 

As  well  by  feeding  back  the  anode  offgas  to  the  burner  and 
the  optimisation  of  the  reformer  insulation  noticeable  efficiency 
increases  can  be  achieved.  Therefore  with  the  tool  it  was  possible 
to  predict  an  average  absolute  efficiency  increase  of  5%  electrical 
and  3%  thermal  for  the  anode  offgas  recirculation  as  well  a  total 
efficiency  increase  of  2.45%  using  a  three  times  better  reformer 
insulation. 

Aside  these  directly  the  existing  system  affecting  achievements 
alternative  system  concepts  can  lead  to  better  total  system  efficien¬ 
cies.  For  the  2  kWel  PEMFC  CHP  system  fuel  treatment  an  alternative 
concept  using  a  selective  methanation  instead  of  a  selective  oxi¬ 
dation  has  been  simulated.  This  concept  allows  a  more  compact 
system  structure,  a  reduction  or  even  removal  of  air  treatment 
and  dosing  components  and  a  simplification  of  the  thermal  system 
concept.  Together  with  the  above-mentioned  direct  system  opti- 
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misation  strategies  this  leads  to  possible  system  net  efficiencies  of 
22%  electrical  and  35.5%  thermal  compared  to  actual  7.7%  electri¬ 
cal  and  21%  thermal.  The  predicted  values  are  in  good  agreement 
with  efficiencies  of  comparable,  real  pre-commercial  systems  and 
prototypes. 

The  developed  modular  model-aided  system  analysis  and  opti¬ 
misation  environment  presented  in  this  work  enables  an  easy 
and  fast  optimisation  of  prototype  systems.  The  tool  also  allows  a 
detailed  characterisation  and  analysis  of  systems  and  components 
and  gives  us  the  possibility  for  flexible  and  dynamic  simulations  in 
the  fast  developing  process  of  new  systems. 
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